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Abstract 

We develop in this paper a new framework for discrete calculus of variations 
when the actions have densities involving an arbitrary discretization operator. 
We deduce the discrete Euler-Lagrange equations for piecewise continuous criti- 
cal points of sampled actions. Then we characterize the discretization operators 
such that, for all quadratic lagrangian, the discrete Euler-Lagrange equations 
converge to the classical ones. 
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1. Introduction 

Discrete variational problems and discrete mechanics constitute active fields 
of research and have been thoroughly studied by many authors among which 
we may cite J. Marsden and J. Cresson, see for instance 0, 0, [H, IE 0] ■ A first 
approach has its roots in the work of J. A. Cadzow, developed in the 1970s, and 
consists in replacing dynamical variables and paths by discrete variables and 
paths which are finite linear combinations of indicator functions. Based on this 
discretization, the action integral, the lagrangian and the energy are approxi- 
mated on a time slice [ke, (k + l)e], for some time delay e, by discrete analogs 
which are finite sums depending of a finite number of variables in M. d . The 
presentation of this theory, together with an history and large bibliography, is 
developed in Q . Another different approach consists in replacing the derivative 
x(i) of the dynamical variable x(t) with a three terms scale derivative D £i gx(t). 
In this way, a dynamical variable x(t) remains even after the sampling process 
a function O e Qx(t) of "continuous" time. The principle of least action may 
be extended to the case of non-differentiable dynamical variables of Holderian 
regularity C 1 / 2 . We refer to [l, 2] for further details among which the extension 
of Noether's Theorem and also an interesting informal discussion about fractal 
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physics. While the first approach is widely used in numerical analysis of varia- 
tional dynamical systems, the second one is, roughly speaking, more concerned 
with foundational discussion of microphysics. 

In this paper we highlight and use the equations of motion 

V x £ + D_ £ Vi£ = 0, (1) 

for all lagrangian C and general discretization operator D s defined as 

JV 

□ £ x(i) = ctx(t + ie)x-t(t). (2) 
e=-N 

In this formula, the symbols X-i{t) denote characteristic functions of some 
intervals, and prevent t+£e from belonging to an interval in which x is undefined. 

A first motivation in this paper is to justify the choice of D £ .q in [H, 0- 
To do this, we use throughout the previous generalization @ and distinguish 
among the operators D e those having some specific features. A second intend is 
to understand the similarities and differences between the equations of motions 
for discretized dynamical variables and the classical Euler-Lagrange equations 
C' x — d/dtC 1 ^ = 0. A third one is to generalize some results of P,0 to operators 
□ e not satisfying some Leibniz formula and to extremal curves x(t) which are 
not C 1 / 2 but rather of a weaker specific regularity C pw . 

In Section 2, we first fix the notation used in the paper and study two 
specific sets On,s and On,s of operators ©. In Section 3, we describe a class of 
operators d e satisfying an appropriate extension of Leibniz formula. In Section 
4, we get the necessary first order condition for finding a minimizer of a discrete 
action Ai«sc(x). We prove the equation (JT]), which depends only on D s and is 
reminiscent of Cresson's result. In Section 5, we first introduce the classical and 
discrete models for the quadratic lagrangians and emphasize on their similarity. 
Next, the subset of operators O s for which discrete Euler Lagrange equations 
have oscillatory solutions is examined. In Section 6, we state the main result of 
the paper which characterizes the set of operators D e for which the convergence 
of ([1]) holds, for every quadratic lagrangian. 

In this paper, C.E.L./D.E.L. stand for classical/discrete Euler-Lagrange 
equations. 

2. Notation and framework 

First, let us fix some notation. Let [a, b] be some interval of time and a time 
delay e > be fixed throughout. We denote by d the "physical" dimension and 
by N the number of samples in C d . We use the notation i for We denote 

by xe{t) the characteristic function of the interval [max(a, a+te), min(£>, b + £e)], 
for all integer I. 

Many interesting ways to generate discrete actions arise in the following way. 
The derivative x(t) occuring in lagrangians is replaced with a linear combination 



2 



of discretized values of x(t) 



□^x(t) = - x+ i(i)-x(t - e) + — x(i) + x_i(t)-x(i + e) (3) 

£ £ £ 

where r,s € C*. In 0, 0], the authors choose the values r — (1 — i)/2 and 
s = (1 + i)/2 but do not work with characteristic functions. Let O e ,Q be defined 

up to now by De . Obviously, ([3]) is a particular case of ([2]) with N = 1, 

s = — c_i£, r = ci£ and Co = — (c_i + ci). We generalize this process by using 
linear combinations of 2N + 1 terms of the shape cix.{t + ie)x-i{t) where q € C 
are fixed. 

Definition 2.1. Given a and (3 in C d , the affine space C pw (d, N, a, (3) is the set 
of the functions x : [a, b] — > C d satisfying x(a) = a, x(o) = (3 and continuous 
on each interval [a + te, a + (£ + l)e] l~l [a, 6] /or integer t. 

This space is endowed with the topology of uniform convergence, its tangent 
space is everywhere the Banach space C pw (d, N, 0, 0). 

Definition 2.2. For any x g C pw (d,N,a, 0), we denote by 5(x) the row vector 
valued function 

S(x)(t) = (an(t - A^e)Xiv(t) • • • a;i(t + te)x-t{t) ■ ■ ■ x d {t + Ne) X -fr(t)). (4) 

In S(x) the ordering of variables is lexicographic, first on £ and next on 
j. The function S(x) lies in a product of d(2N + 1) affine spaces modeled on 
C pw (l, N, 0, 0) and, in this way, S is an injective continuous linear mapping. 

Definition 2.3. If(ct) £ C 2N+1 , the generalized scale derivative is the continuous 
linear endomorphism of C pw (d, N, a, j3) defined by (0) for all t € [a, 6]. 

If d = 1, for all x e C pM ,(l, iV, a, /3), we have 

□ e x = (c_jv • ■ • c . . . cat) T (5(x)). 

When d > 2, we have a similar relationship between D e x and S'(x), involving a 
d x 0^(2^ + 1) banded matrix 



CL-X : 



/ C-AT ... C ... Cat ^ 

c-n ... c ... Cat 



\ C-N ...C ...Cjv/ 



T (5(x)). 



Note that D s is a well defined continuous endomorphism of the Banach space 
C pw (d,N,a,/3). 

The convergence of D e x to the ordinary derivative x is connected to the following 
definitions. 

Definition 2.4. The vector space On,s is the set of operators D e of the shape (0) 
with coefficients ci — ji/e and where 7^ G C do not depend on e. Algebraically, 

o Ni£ ~c 2N+1 . 
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Definition 2.5. The affine space On,c C On,s contains operators \3 e for which 
□ e l = and D e t = 1 when t lies in the safety interval Is = [a + 2Ne, b — 2Ne] . 



We shall see, as a consequence of Theorem 16.11 that if O e 6 On.e, then 
□ e x(£) tends to x(t) locally uniformly in }a, b[, for all x £ C 2 ([a, 6], C d ). 

Remark 2.1. The discrete Euler forward and backward difference operators are 
respectively defined by A+ = □J, 1,0 ' and A~ = d! ' 1 '. Their mean is the sym- 

metric difference operator Sj s = D/' 2 • These three operators and O e ,Q are 
elements of 0\. e . Notice that the operator O e ,g is also related to the one- 
dimensional version of the operator used by Kime [J] when she solves numerically 
the Schrocdingcr equation via Forward-Difference and Backward-Difference ap- 
proximations at various steps in time. 

Remark 2.2. The operator x — > A+(A+x(i — e)) = — □[ 1,_1 'x(i) approximates 
well the second derivative but does not lie in 0\_ E . 



3. Leibniz formulas for d e operators in 0± tS 

In order to deduce his version of D.E.L. with D^g, Cresson in [1| found an 
analog of the classical Leibniz formula. When such a formula exists, a principle 
of discrete virtual works may be stated. In this section, we generalize Cresson's 
identity to the family of operators 

Theorem 3.1. Let O e € Ci, e of the shape f3|) with r, s € C* and s/r ^ ffi. then, 
for all piecewise continuous functions f,g : [a,b] — > C d , we get the generalized 
Leibniz formula 

□l^i (f ■ g)(t) = f(t) • nM g (t) + g (t) • □r^]f(t)+ 

e(rs 2 — r 2 s) 



;(*)- 



(rs — rs) 2 e e (rs — rs) 2 

(4 M fW • □P ] gW + np ] fW ■ □i r ' s] g(0) • (5) 



Proof. The formula ([5]) is C-bilinear w.r.t. f and g. Having in mind the prop- 
erties of the inner product and the fact that D e acts component-wise, we may 
suppose without loss of generality that d — 1 and f , g G R. We slightly gen- 
eralize the proof of Theorem 2.1 of [l( which attempts to give a formula such 

as 

W(fg) = W(f)g + fW(g)+ 
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dxW{f)W{g) + d 2 W(f)W(g) + d 3 W(f)W(g) + <kW(f)W(g). (6) 
Here W and W are two operators in Oi. £ , f and g are arbitrary in C pw (l, 1, a, /3) 
and ^1,^2,^3,^4 are four complex numbers. Now, we choose = □]T ,S ' and 
W = ,s ' for some complex numbers r, s, r', s' such that rs' — sr' ^ 0. We 
have obviously W = rA+ + sA~ and = r'Af + s'A~. But we have also the 
well-known formulas for Aj~ and Ar 



A°(fg) = AUf)9 + fK{g) + veA°(f)A°(g), 



(7) 



where a — ±1. Substituting the four previous formulas in ([5]), the identity ([5]) 
holds if and only if the coefficients d\,di-> d$ and di satisfy 



/r 2 



rs rs 



V 



sr 
ss' 



o'2 



\ 






(re\ 













d 3 







/ 


{d 4 J 




{-ssj 



(8) 



whose determinant is equal to — (rs'— sr') 4 ^ 0. We replace in (jSJ) the coefficients 
di by their explicit values 

di = |(rs' 2 - sr' 2 ), d 2 = d 3 = - s'), d 4 = ^(s - r) (9) 

where <5 = (rs' — sr') 2 . Since s/r ^ K, we can choose r' = r and s' = s and we 
get easily the formula ([5]). □ 

As an example, we get —di = d 2 = d 3 = d± = — \ie for the operator O e ,Q — 
□e 2 ' 2 ' chosen in [H, Q, so that for all piecewise continuous /, g : [a,b] —> K 

□ e ,q(/ 5 ) = □ £ , S (/)<7 + /a,s(5) + \ie[U e>Q {f)n e>Q {g) - U e , Q {f) B £ , s ( 5 ) - 
B E ,e(/)n e ,e(0)-B e>a (/) B E>S (<?)], 

where B e stands for the complex conjugate operator of D e . (We have corrected 
here the corresponding formula in [l|, df.) 



4. Critical points of discrete actions 

According to the previous notation, the discrete actions and lagrangians are 
related to classical ones as follows. 

Definition 4.1. If C(t, x(t),x(t)) is a lagrangian depending on 2d + 1 variables 
then we set L(t, S(x)(t)) = C(t, x(i), O e x(t)). Moreover, 



and 



A cont (x)= / C(t,x{t),i(t))dt (10) 

J a 

rb rb 

A disc (x)= L(t,S(x)(t))dt= £(t,x(t),a £ x(t))dt (11) 

J a J a 

are the respective classical and discrete actions associated to these lagrangians, 
defined respectively on C^fa, b], C d ) and C pw (d, N, a, f3). 
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If the terms Xj(t + £e)x~e{t) occuring in S(x)(t) (see (|3])) are replaced with 
new variables then L(t, £i,-jv, . . . , £d,iv) is nothing but £(t, x(i), D e x(i)) 
and thus depends on d(2N + 1) + 1 indeterminates. 

A fundamental problem is to minimize the action (jlOp under Dirichlet bound- 
ary conditions when obviously, every unknown and parameter has to be real. 
Note that in order to deal with optima instead of critical points of the action 
(jlip . we have to handle real valued functions and parameters. 

Theorem 4.1. Let x € C pw (d, N, a, /?) be a critical point of Then x satisfies 

the following functional equation 

- 3J 

Vje{l,...,d}, V (t-£e 1 S(x)(t-ee))xi(t)=0 (12) 

i=-N °W 

and the equation Uty) may be returned under an intrinsic form as (QP, i.e. 

||(t, x(t), □ e x(t)) + ^-(t, x(i), ax(t)) = 0, (13) 

/or a// j G {1, . . . 

Proof. We have for all h e C pti! ((i, N, 0, 0) 

Afec(x + rjh) - A dlsc (x) = f (L(t, S(x + nh)(t)) - L(t, S(x)(t)))dt 

J a 

.b d N „ 

= 1 EE ^—{t,S^){t))h j {t + is)x-t{t)dt + 0{v 2 ). 

Therefore using Chasles relation and setting t = t — is in the previous equality, 
we find the Gateaux derivative 



d n r b+i s g L 



D^ sc ( X )(h) = EE / 

7=i e=-N , ' a 



fb+Nc d N „ 

/ E^'M E ^-(r-fe,S(x)(r-&))x<(r)dr = 0, 



(t - ie, 5(x)(r - te))hj{r)dT 



which gives (1121) by using Paul Dubois- Reymond lemma in the multidimensional 
case. By definition of L(t, 5(x)(t)), we have for all £ ^ 0, 

<>L (i,5(x)(t)) - !^(t,x(t),n e x(*))c/X-< 



3^ ctej 
and for € = 0, 

— - t, S x i = — (i, x(t), ax(t))c + ^— (t,x(t), n e x (t . 
c^o oxj axj 

Thus, (p~2]) is equivalent to 



G 



N dC dC 

c «-^(* - fe > - D ^ x ( < - fe )) + ^t(*> x W' Q ex(t)) = 0. 



i=-N 



dC 

Since the first term is equal to D_ E - — (i, x(t), D E x(t)), the result is proved. □ 

Remark 4.1. When D e is not of the shape we do not have and integration 
by parts may not be performed. We use instead simple changes of variables. 
Note also that second order derivatives of the lagrangian occuring in C.E.L. are 
replaced with time delayed first order ones. 

Remark 4.2. In [jj, Cresson deals with the case N = 1 and the operator Si q = 

□e 2 ' 2 . For any function /(x, e) he defines the e-dominant part [f] e with a 
limiting process. He proved that if 



lim DAdisc(*) = 0, 

E— >0 



then 



>6 \ dxj e ' q dx 3 



^(—-□.,^1=0 (14) 



which is an alternative form of (1131) . Our equation (jl3[) is however different 
from (fT4]) at least in three respects. First, the characteristic functions X-e(t) of 
the various intervals appear in the sampling process ([2| as well as in the action 
pT|) . The second one is that (TlUl) does not depend on the coefficients nor on 
the length of formula @. The last one is the use in (IT5|) of D_ e instead of 
— n s . But we have [□ e) g/(t)] E = [—□_£. g /(t)] e for all real valued function / 
and t G [a + e, b - e}. Indeed, [S(D e ,s/j] e = and 9t(n_ e>q f) = -R(D e>s /) 
since every characteristic function in De.g is equal to 1. 



5. Quadratic lagrangians in discrete and classical settings 

In this section we deal with a system of d ordinary differential equations of 
the second order arising from the following lagrangian 

C(t, x,x) = - xPx + - T xQx + VRx + T Jix + T J 2 x + J 3 , (15) 

where P(t),Q(t),R(t) G C dxd , Ji(t),J 2 (t) G C d and J 3 (t) is a scalar function. 
Many physical systems might be modelized by such lagrangians, in electromag- 
netism, quantum mechanics, material science, regulators models and so on. 

A convenient setup that we assume from now on is that the coefficients in 
(fT5|) are real and smooth, and that for all t G [a, b], P(t) and Q(t) are symmetric 
and R(t) is skew-symmetric (the symmetric part of R(t) gives rise to a null 
lagrangian term in (|T5|) ). 
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Theorem 5.1. Let C be a quadratic lagrangian such that T P = P, T Q = Q, 
T R = —R and L associated to L as in ill)) . The Euler- Lagrange equation 
associated to U5\) can be written as 

- Px + (-P + 2F)x + (R + Q)x - j x + J 2 = 0. (16) 

The equation M6\) may be discretized a posteriori to give 

- ra £ (D E x) + (-P + 2R)D £ x +(R + Q)x - J x + Ji = 0. (17) 
IfxE C pw (d,N,a,j3) is a critical point of the action Ml)) , then it must satisfy 

□_ e (PD E x) - CL E (Px) + PD E x + Qx+ □_ E Ji + J 2 = 0. (18) 



Proof. First, (|TB)) is straightforward, since we get 

= Qx + i?x + J 2 and = Fx - i?x + J, . 

ax ax 

Next, (|T7| is obtained by discretizing the derivatives in (fT6|) . i.e. by replac- 
ing x,x by D E x, □ e (D e x) respectively and the result holds. At last, (fT8|) is a 
consequence of (|13jl . Indeed, using the previous derivatives, (|13|) gives 

□_ e (PD e x - Px + Ji) + Qx + RD E x + J 2 = 0, 

which ends the proof. □ 

Remark 5.1. If P, Q, R are independent on time then, for t € Is, both equations 
(fT5|) and (|17[) are equivalent. 

Remark 5.2. For the so-called linear quadratic control problem, Guibout and 
Scheeres Q give Euler-Lagrange differential equations having a similar structure 
than (HU). They suppose that the relation x(t) = M x (i)x(i) + M 2 (i)u(i) holds 
between the dynamical variable x and the control variable u and prove that 
each critical point (x, p) satisfies 

x{t) = (Mi - \AhP- 1 T R)x{t) - M 2 P- X T M 2 p{t), 
p(t) = (iP T P- x T Af 2 - T M!)p(i) + (IPP- 1 T P - Q)x(t), 

where p stands for a Lagrange multiplier of the constraint. 

Let us investigate now the harmonic oscillator. We have in this case d = 1, 
R = 0, Ji = J 2 = and we set moreover P(t) = p, Q(t) — q with pq < 0. 
Solutions of C.E.L. —px + qx = are periodic. Let us study the periodicity of 
solutions of D.E.L. for e small enough. If t lies in Is, ([21]) may be simplified 
into 

pc_ici[x(< - 2e) + x(i + 2e)] + pc (c_i + ci)[x(t - e) + x(t + e)] + 

(g + p(c 2 _ 1 + C 2 + c ?))x(t) -0. (19) 

The result below presents an additional characterization of "suitable" O e . 
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Proposition 5.1. Let \3 S £ Oi, £ . The two following properties are equivalent. 

(a) For all p and q with pq < and for all e small enough, the roots of the 
characteristic polynomial of H9\) are of modulus 1. 

(b) For some k G R, we have D £ = D^ 1 + ikn£'~ l] . 

Proof. The characteristic polynomial of the recurrence (| 19[) is symmetric that 
is to say D(X) = X 4 D(1/X). Setting fj, = X + 1/A, we get a quadratic E(/j) such 
that 

E{p) = e 2 D(X) = 7i7_i/i 2 + 70(71 + 7-i)P + (To + 7? + 7-i - 2717-1 + -e 2 )- 

P 

(a) (6). We look for the parameters 7-1,70,71 for which D(X) has roots on 
the unit circle, for e small enough and for all p, q, pq < 0. This amounts to say 
that E(fi) has only real roots in [—2,2] for all p,q, pq < 0, provided e is small 
enough. Let us generalize a little bit. Let a, /3, 7 be three complex numbers, 
then the solutions of the quadratic ay 2 + [3y + j + e — are in [—2,2] ct for 
all e small enough if and only if we have a, j3, 7 G R and 

< 4|a|, |7| < 4|a|, 4a 7 < /3 2 , 8|/3| < 16|a| + A-fsgn(a), 

as shows explicit computations using the usual solution of a quadratic. Since 
□ e € 0i, 6 , we have O e l = 7 _i +70+71 =0 and D e t = 71 — 7-1 = 1 inside 
Zs if and only if 71 = r, 7 _i = r — 1 and 70 = — 2r + 1 for some r € C. So, 
a = r(r — 1), /3 = (2r — l) 2 and 7 = 4r 2 — 4r + 2 are real. Therefore we obtain 
K(r) — 1/2 and easy computations lead to (6). 

(6) (a). If (b) holds then setting r = \ + jfe, we get 71 = r, 7 _i = 1 — 1 and 
70 = — 2?* + 1. Direct computations show that the solutions of 

E(p) = r(r - l)^ 2 - (2r - l) 2 // + (4r 2 - 4r + 2) - uj 2 e 2 , (20) 

where tu 2 = -q/p, are in [-2,2] c M for all e < 1/(M%/1 + 4/c 2 ). □ 

Remark 5.3. We recover D^g and De.s by setting fc = —1/2 and fc = in (6) 
respectively. 

Remark 5.4. If D e lies in Oi, e , the property (a) is equivalent to the following 
inequalities for 7-1,70,71 G C 

4|7i7-il - |7o(7i +7-i) I > 0, 4| 7 i7_i| - I7 2 +7? +7^i - 2 7 i7-i| > 0, 
7o(7i +7-i) 2 - 4717-1(7^ +7 2 + j 2 _ 1 - 2 7 i7_i) > and 
I6I717-1I +4 7 2 + 4 7l 2 + 47^! -8717-1 -8|7o(7i +7-i)l > 0. 

So there exists operators D e ^ 0\^ e satisfying (a) but not (b). 
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6. Convergence of functional equations D.E.L. to differential equa- 
tions C.E.L. 

In this section, we address the problem of convergence of D.E.L. to C.E.L. 
in the case of quadratic lagrangians. Roughly, this convergence property char- 
acterizes Ojv, £ - In the remainder, the matrices P, Q, R in (|15p are dependent on 
time, P, Q are symmetric and R is skew-symmetric. To begin with, let us give 
the following 

Lemma 6.1. The left hand side of D.E.L. \18\) is equal to 

0(x)(i)= c £+j c jXj (t)x-e{t)P(t-jeMt + £e) + Q(t)xL{t)+ 

~2N < t < 2N 
-JV < j < N 
JV \t + j\ < N 

J2x-e(t)(c e R{t) - c- t R(t + fe))x(t + te) + n_ e Ji{t) + J 2 (t) (21) 

e=-N 

for all t S [a, b]. 

Proof. The proof is straightforward by using several times the formula © for 
□ e and □-£, and showing that (fT8]l may be returned as 0(x)(i) = for all 
t e [a,b]. □ 



Definition 6.1. We say that D.E.L. U8]) converges to C.E.L. U6]) as e tends to 
if, for all quadratic lagrangian C as in U5\) and all x € C 2 ([a, b], K d ), 

lim 0(x)(t) = -Px + (-P + 2P)x + (R + Q)x - j a + J 2 

locally uniformly in ]a, 6[, m £/ie norm of C°°([a + 5,b — 8]) for all 8 > 0. 

The remainder of this section is devoted to the proof of the main result 
below. 

Theorem 6.1. Let d e € On.c- The following six properties are equivalent. 

(a) D.E.L. U8\) converges to C.E.L. \16\) as e tends to 0. 

(b) For all x G C 2 ([a, b],M. d ), Una □gx(t) = x(t) locally uniformly in ]a,b[. 

(c) For all x £ C 2 ([a, b], M. d ), lim D_ e x(t) = — x(t) locally uniformly in }a,b[. 

s— >o 

(c?J TTie functions t — > D e l and i — > D e t converge respectively to and 1 
locally uniformly in ]a, 6[. 

f/y There exists complex numbers ki, . . . , fejv-i suc/i £/ia£ 

2V-1 

□ £ x(t) = □I 1 ' 0] x(t) + ^ kt+N^'-^xit - &). (22) 
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Proof. We will prove that (a) =► (c) =4> (d) (e) ^> (/), (d) =► (6) (c) =► (a), 
(a) (c). By assumption, the function 9(0) = D-gJi + J2 must tend to the 
constant term in (fTo| that is — J\ + J2, for all J\ and J2, and the result holds. 

(c) =>■ (d). We choose x as a linear function t — > Xii + Xo which checks the 
property in (c). As D_ e acts component- wise, D-^l and [ZL e f tend to as s 
tends to 0. But since 

1 N 1 N 

□ e l = - J! 7a-/(*)andn e t = tD e l + - t(leX-t(t)-l-eXi(t)), (23) 

for all £ 7^ 0, we have inside 1$ the relations O e l = — D_ e l and D e t = — D_ e t + 
t(D e l + n_ e l). This shows that (d) holds. 

(d) (e). For all S > and £ small enough, if t S [a + S, b — S] then t £ Zg, and 
we see that the assumptions in (d) are equivalent to (I23p and in turn to the two 
linear equations 

jv 1 N 

J2 It = and - J] £( 7£ - 7 _<) = 1. (24) 

£=-iV l=-N 

If n e <E Ojv,e satisfies these two equations then D £ l — > and D £ t — > 1 uniformly 
locally in ]a, b[, that is O e satisfies the statements in (d) and conversely. Thus, 
we have proved that D £ £ On,c if and only if (d) holds. 

(e) (/). The two equations (|24|) being independent, codim(0jv :£ ) = 5. If 
On,c is the set of operators defined by (f22j) . then codim(Ojv, e ) = 2- We easily 
see that each operator D E £ On :£ checks the two equations in (d), that is 
On,e C On.e and lastly, On,e = On, e - 

(d) (b). We use Taylor-Lagrange formula to obtain an expansion of D e x as 



JV 



t+ts 

x(i) + fex(t) + / (t + £s- s)x(s)ds 



□ e x(t)= £ c«_*(t) 

= (□ e l)x(t) + (D e < - tD e l)x(t) + / G(s, i)x(s)ds. 

J a 

By using X[o,fc] = X[o,fc] if ^ > and X[o,&] = X[£e,o] if ^ < 0, the previous kernel 
G(s, t) is equal to 

N 

G(s, t) = X o(s)xo(t) ^(t + ie- s ) c eX-e(t)X[o,ee] - t). (25) 
t=-N 

Let 5 > 0, V x = C°°{[a,b]), V 2 = C°°([a + S,b - 5}) and V 3 = C°°{[a,b] x [a + 
5, b— 5]). Then, for all function x £ C 2 ([a, b]), if £ < <5/iV, the norm ||D E x — x\\v 2 
is bounded by 

||n E l||v 2 (||x|k + max(|a|, |6|)||x|| Vl ) + ||D e i - l|k||x|k + (b - a)\\G\\ Va \\i\\ Vl . 
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The first and the second terms converge to as £ tends to help to (d). In 
order to prove that the last term tends also to as £ tends to 0, we note that 

N 

G(s,t) = (t-ap e l+ ( fc X[o,fe](5-t) + (t-s)(x[o,fc](^0-l))c«-£W- 
t=-N 

Obviously, first and second terms tend to as £ tends to 0. To deal with the 
third term, we distinguish two cases. If t < s < t + fee for some k € [— N, N], 
the term is bounded by Ne^2f = _ N \ce\ while if \s — 1\ > Ne, it is bounded by 
\s — t\.\\D £ l\\v 2 , and thus, by the sum of these two upper bounds. Collecting 
altogether, we find 

N N 

\\G\\ Vs <2max(|6- a |,2|a|,2|6|)||n e l||v a +JVe £ \c e \ + e £ \lc t \. 

t=—N t=-N 

Therefore, ||(?||v 3 tends to and thus, also ||d e x — x||yi as e tends to and the 
result holds. 

(b) =$> (c). The two operators O e and n_ e are linked by the following property. 
If xi <E C pw (d, N, a, (5) then the function X2 defined by x 2 (t) = Xi(a + b — t) is 
in C pw (d, N, /?, a) and satisfies 

(□- e xi)(*) - -(n e x 2 )(a + b-t). 

Indeed, this formula comes from ((2} and the fact that x-j( s ) — Xj( a + b — s) 
for all j. Since ^(t) — — x'i(a + b — t), the convergence of D e X2 to X2 for 
all X2 £ C pw (d,N,/3,a) implies the convergence of D_ e xi to — xi for all xi G 
C pw (d,N,a,/3). 

(c) =£> (a). We use again the formula (|2"TT) and we deal with its non-constant 
terms. When we develop x(£ + £e) help to the Taylor-Mac Laurin formula of 
order 2, we get 

9(x)(t) = -P £ (i)x(t)+2 J R e (t)x(i)+Q E (t)x(t)+n_ e J 1 (i)+J 2 (t)+9(x)(t), (26) 

where notation is explained hereafter. First, let us note that 0(x)(i) is a com- 
bination of two kinds of terms M(x(i + — x(t)), for some #g €]0, 1[ and 
€ e {-2iV, . . . , 2iV}, M being either equal to P(t - j'e) or R(t + le). When 
t € Ig, we see that those terms tend to as £ tends to 0. Next, some lengthy 
computations show that 

p £ = -±p[a- £ n e t 2 -2ta- e n e t + t 2 a^ £ a £ i] + s P , 

Q e = Q + pn_ s n s i + p(n e i + n_ e i)-p[(n e i)(n s t-tn e i)]+ 
P[(D e i)(a e t 2 - 2tn e t + t 2 n s i)] + R(-n_ e t + <o_ e i) + s Q , 

i? e = ip(D_ e n e t-tn_ e n e i) + ii?(n e t-n_ e t + tn- e i-tn e i) 
-§p[(n B t - tn E i) 2 - (□ e i)(n e i 2 - 2ta E t + t 2 n e i)] + 5 fl . 
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Sp, Sq and Sr are three explicit matricial combinations of P, P, P, R, R, R 
which may be roughly upper bounded as follows 

INU~(M]) < l 2 N(e + x^) (E \\P {e) \\c-aa M ) + E P W IU~(M])j • 

Here, 7 = max(l, max< |7^|) and N = (2N + l)(4iV + 1)N 5 are two fixed numbers 
w.r.t. e. As a consequence, for all small S > 0, XxjW an d accordingly <Sj(£) 
tend to uniformly in [a + 5, b — 5} as e tends to 0. Finally, since (c) holds, 
inspection of each coefficient in the previous formulas shows that P s (t), Qe{t), 
Re(t) and Q(x)(i) tend respectively to P(t), Q(t) + R(t), R(t) - \P{t) and 
as e tends to 0. The convergence of D.E.L. is thus ensured for all lagrangian £, 
for all function x and for all t € [a + 6, b — 8] , which ends the proof. □ 

Remark 6.1. We note that inside the safety interval Ts, for each operator 
□ e G O n ,s, the three formulas hold: U- e U e l = (D £ l) 2 , U- e U e t = t(D e l) 2 

and a_ e D e t 2 = t 2 {D £ l) 2 + 2(D e 1) I E _ 2 ( D e*) 2 - More generally, 

\i=-n ) 

there exists polynomial formulas for iterates of D E acting on the polynomials as 
expressions of \3 e t k for fc e N. 

Remark 6.2. As mentioned in Remark 14.21 Cresson defined extremal curves 
of the action as functions x(t) such that [£>^l [ ;isc(x)] £ = 0, that is to say the 
dominant part of the Frechet derivative of the action is zero. Obviously, this 
does not ensure that Adisc is extremal at x. The matter of the convergence of the 
Euler-Lagrange operator restricted to extremal curves of Holderian regularity 
C^([a — e, b + e],M. d ) for some (3 > forms part of the characterization of 
Discrete Euler Lagrange equations. In contrast, we work here with all curves 
x 6 C 2 ([a, b], R d ), not necessarily extremal. 

Remark 6.3. In the terminology of T- convergence, see for instance Q, the pre- 
ceding proof implies that for all t Q e]a, b[, if we denote ev to the evaluation map 
at to, then the composite mapping x — > ev to o 0(x) is T-convergent to 

x -> ev to o (-Fx + (-P + 2R)± +(R + Q)x - A + J 2 ). 

Note that the authors in [7| proved the T-convergence of actions which are 
quadratic w.r.t. x and defined on some spaces of piecewise affine maps. 
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